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Abstract 

The relativistic nuclear recoil corrections to the energy levels of 
low-laying states of hydrogen-like and high Z lithium-like atoms in all 
orders in aZ are calculated. The calculations are carried out using the 
B-spline method for the Dirac equation. For low Z the results of the 
calculation are in good agreement with the aZ -expansion results. It is 
found that the nuclear recoil contribution, additional to the Salpeter's 
one, to the Lamb shift (n = 2) of hydrogen is —1.32(6) kHz. The 
total nuclear recoil correction to the energy of the (ls) 2 2pi — (ls) 2 2s 
transition in lithium- like uranium constitutes — 0.07 eV and is largely 
made up of QED contributions. 

PACS numberfs): 12.20.Ds, 31.30.Jv 
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1 Introduction 



As is known, in the non-relativistic approximation the nuclear recoil cor- 
rection for a hydrogen- like atom can be taken into account by using the 
reduced mass /x = ™+ M - The relativistic corrections of order («Z) 4 ^mc 2 
can be found by employing the Breit equation [1]. A theory of the nuclear 
recoil effect in higher orders in aZ must be constructed in the framework of 
quantum electrodynamics (QED) on the basis of an exact relativistic equa- 
tion for hydrogen-like atom. Such an equation was proposed by Bethe and 
Salpeter [2] immediately after creation of QED. On the basis of this equa- 
tion the nuclear recoil corrections were calculated in [3] up to terms of order 
(«Z) 5 |mc 2 . It was shown in this work that the nuclear recoil effect in the 
case of a complex nucleus is calculated in a good approximation by assuming 
the nucleus is the Dirac particle with the charge \e\Z and the mass M. Sub- 
sequently these corrections were recalculated by a number of authors [4-6]. 
Calculations of the nuclear recoil corrections of the next order in aZ were 
considered in [7-11]. 

In the theory of high-Z one-electron ions the parameter aZ can no longer 
be considered small. For this reason calculations of the nuclear recoil correc- 
tions for such systems must be carried out without expansion in aZ. In con- 
trast to other QED effects in the region of strongly bound states (aZ ~ 1), 
the calculation of the nuclear recoil effect at high Z demands using QED 
outside the external field approximation. ( Calculations of QED effects in 
hydrogen, positronium, and muonium correspond to the case of weakly bound 
states (aZ <C 1).) In this connection a non-trivial problem of derivation of 
closed expressions for the nuclear recoil corrections in all orders in aZ arises. 
This problem was first discussed in [12,13]. The work [12] was based on 
the Bethe-Salpeter equation. This approach encountered serious technical 
difficulties, associated with summation of a complete sequence of irreducible 
diagrams. These difficulties were partly overcome only in the lowest orders 
in aZ. Complete aZ-dependence expressions were not found in this way. In 
[13] a general case of a relativistic few-electron atom was considered. An ef- 
ficient method for summing of the Feynman diagrams in the zeroth and first 
orders in based on an expansion of the nuclear propagator, was proposed 
in this paper. However, because the procedure of the derivation of the nu- 
clear recoil corrections was not rigorously formulated, the method considered 
there gave several ambiguities in the expressions for the nuclear recoil cor- 
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rections. In addition, certain errors were made in derivation of the formulas 
for the contributions with one and two transverse photons. As result, only 
a part of the expressions for the relativistic nuclear recoil corrections was 
found in this work. The complete expressions for the nuclear recoil correc- 
tions for hydrogen-like atoms were obtained in [14] (the overall sign of the 
two-transverse-photons contribution was corrected in [15,16]). The paper 
[14] was based on a version of the quasipotential approach that immediately 
gives the Dirac equation in the limit of infinite nuclear mass [17,18,5]. ( 
The quasipotential approach was first introduced in quantum field theory by 
Logunov and Tavkhelidze [19] and was subsequently developed by many au- 
thors (see, e.g., [20]). This approach is absolutely rigorous and, in contrast to 
the Bethe-Salpeter equation, allows one to exclude the relative time (energy) 
in the wavefunction from the very beginning. The quasipotential equation 
can be represented in the evidently covariant form [20,17].) The relevant 
quasipotential equation in the center-of-mass system is (the relativistic units 
h — c — 1 are used) 

(E - Vp 2 + M 2 - ap - /3rr#(p) = J V(E, p, q)^(q)rfq , (1) 

where ex., (3 are the Dirac matrices acting on the electron variables. The 
quasipotential V(E) can be constructed by various methods [17,19,20]. One 
of the methods consists in using the relativistic scattering amplitude with one 
particle (nucleus) on mass shell [17,18,21]. In this method the quasipotential 
V(E) may be defined by the Lippman-Schwinger equation 



where 



V = t(1+Ft)-\ (2) 

F = [E - (yp 2 + M 2 + ap + /3m) (1 - ^O)]" 1 , (3) 
r(£,p,q) = -2mPn(-p)T(p 1 ,p 2 ;q u q 2 )u(-q) , (4) 

Pi = -P2 = p , q i = -q 2 = q , 

E- yV + M 2 , pl = vV + M 2 , 



Pi 



q» = E-^q 2 + M 2 , q^^cf + M 2 , 

Pi, qi are the electron variables, p 2 , q2 are the nucleus variables; T is the off- 
mass-shell relativistic scattering amplitude; w(q) is the wavefunction of the 
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free nucleus with the positive energy normalized by the condition w(q)u(q) = 
1 . In [14] the quasipotential V(E) was constructed in the zeroth and first 
orders in jj. So, the closed expressions for the nuclear recoil corrections in 
the first order in ^ and in all orders in aZ were obtained. The most detailed 
derivation was published in [22]. In [16] these results were generalized to the 
case of high Z few-electron atoms. For that a more general method was 
developed. In the second section of the present paper we briefly formulate 
the results of [16]. In the third section the calculation of the nuclear recoil 
corrections for hydrogen-like atoms is considered. In the fourth section the 
corrections for high Z lithium-like atoms are calculated. 

2 Basic formulas 

We consider the system of Dirac particles: a nucleus with mass M and N 
electrons with mass m. Following to ideas of the quasipotential approach we 
introduce in the center-of-mass system the two-time Green function with the 
nucleus on the mass shell 

G(t', t, p', Xi, ...X^, p', X 1; ...Xjv) 

= <p', \\Tip{t' , xi) • • • ij(t>, ^' N )^(t, XJV ) • • • ^(t, Xl ) |p, A) , (5) 

where ip(x) is the electron-positron field operator in the Heisenberg repre- 
sentation, T is the time ordered product operator; 

|p, A) = a in (p, A) |0) , |p', A) = a out (p', A) |0) (6) 

are the in and out states of the nucleus; p and A are momentum and polar- 
ization of the nucleus. Here we normalize the operators a in and a out by 

{4»(P> A )> a Up', A')} = {4ut(P, A )> a ut(p', A')} = S xy S(p - p') . (7) 
Let us introduce the Fourier transform of G: 

S(E - E')5(P - P')G(E, p', p[, ...p' N , p, Pl , ...p N ) 

dtdt'dx.1 ■ ■ ■ dx N dx.[ ■ ■ ■ dx.' N exp [i(£'t' — St)] 

N 

x exp [-i ^(p-x^ - piXi)]G(t', t, p, xi, ...x^, p', xi, ..x N ) , (8) 

i=i 



2m N\ (2n) 3N 
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where 

E = S + yp 2 + M 2 — M , E' = S'+ yp' 2 + M 2 - M , 

AT AT 

P = p + 5>, p' = P ' + Ei (Q) 

1=1 1=1 

In the center-of-mass system we have 

N N 

p=-Eft. p' = -Ep: 

i=i i=i 

Let we are interested in the energy of a bound state n of the atom. The 
spectral representation of G(E) gives 

_ $ $ t 

G(E) = " + terms regular &t E — E n , (10) 

where _E n is the bound state energy with the nucleus rest mass subtracted, 
the wavefunction $ n is defined by equation 

(27r)§<J(P)$ n (p,p 1 ,...p JV ) 

1 1 r N 
= ncfrTZ ,m d *i--- d *N exp [- ^ p^] 
Viv! (27r) 2 •/ i=1 

x<p|^(0,x 1 )...^(0,x JV )|n). (11) 

The Green function G(E) is constructed by perturbation theory after tran- 
sition in (5) to the interaction representation. Let the energy level n belong 
to a m-fold degenerate level E^ in the limit M — > oo if the radiative and 
interelectronic interaction corrections are neglected. (The neglect of the in- 
terelectronic interaction in the zeroth approximation is justified for high Z 
few-electron atoms (N <C Z).) The m-dimensional subspace generated by 
the unperturbed eigenstates making up this level we designate as Q. The 
projector on Q is 

m 

Po = E«*«t> (12) 



k=l 



where 



u k = -L= (PI) • • -^ kN {PN) , (13) 



p 
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ipk are solutions of the Dirac equation in the Coulomb field of the nucleus: 
Hipk = tk^k , H = ctp + (3m + V c , 

N 

= E e *i- (14) 
i=i 

Let us introduce the Green function g: 

g = P GP . (15) 

For this Green function, like G(E), we have 

g(E) = + terms regular at E = E n , (16) 

E — E n 

where n = P^ n belongs to the subspace Q. Constructing g{E) by the 
perturbation theory in the interaction representation we get it in the form 
of a series in powers in aZ. However, we are interested in an expansion in 
another parameter, namely, ||. For this reason it is necessary to sum infinite 
sequences of the Feynman diagrams in the zeroth and first orders in -p. We 
designate the contribution of the terms of the zeroth order in ^ by go. In 
[16] it was found 

From the equation (16) and the identity 

g-'g = i (18) 

we obtain for E = E n 

g-\E n )<j> n = Q. (19) 
Or, introducing the quasipotential operator 

v(E) = g, 1 - g- 1 = g, 1 Ag g Q l + g, 1 Ag g l Ag g' 1 + • • • , (20) 
where Ag = g — g , we obtain 

{EW+v(E n ))4> n = E n 4> n . (21) 



6 



It follows the equation for determination of the energy levels 

det{(E - E®)6 ik - Vik(E)} = . (22) 

It should be stressed that equation (22) is absolutely rigorous within QED 
and gives, in principle, the exact energies of the m levels arising from the 
m-fold degenerate level In [16] the quasipotential v ik was found in the 
first order in ^ and in the zeroth order in a (but in all orders in aZ) by 
summing infinite sequences of the Feynman diagrams in the Coulomb gauge. 
For that the expansion of the nuclear propagator from [13] was used. Only 
the following kinds of the diagrams contribute in the considered order: 

• The diagrams with only Coulomb photons. 

• The diagrams with one transverse and arbitrary number of Coulomb 
photons. 

• The diagrams with two transverse and arbitrary number of Coulomb 
photons. 

The contribution from the diagrams with only Coulomb photons is 

(v c ) ik = (vi% k + {vf\ k + (vi in % k , (23) 

TV 2 

= ^ ilkl -.-n...5 iNkN (i s \^\k s ), (24) 

2%i N s r°° 

H 2 ')ik = -j^-J2 5 iik 1 ---n---5 iNkN J^duj5 2 + (uj) 

x(i s \[p s ,v s }G s (uj + e is )[p s ,v s }\k s ) , (25) 
1 

M 



s<s' 

X ' 



Y,(-l) P (P*sPts'\PsPs>\k s k s ,), (26) 



p 



where \i s ) and \k s ) are the one-electron unperturbed states of the Dirac elec- 
tron in the Coulomb field of the nucleus, belonging to the N-electron states 
i and k, respectively; p is the momentum operator, v s = V c (r s ) = — ^ ; 

the symbol n means that the factor 5 iskg is omitted in the product; S + (uj) = 
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j-(u;+iO) -1 , G(uj) = (uj — H^l—iO))^ 1 is the relativistic Coulomb Green func- 
tion. ( Formally, the matrix element in equation (25) at fixed u is infinite, 
due to the strong Coulomb singularity at r = . It means that the integra- 
tion over uo must be carried out on an intermediate stage of the calculation, 
depending on which representation of G is used.) The contribution from 
the diagrams with one transverse and arbitrary number of Coulomb photons 
consists of two terms. The first term depends on the spin of the nucleus and 
coincides with the Fermi-Breit expression for the hyperfine interaction [23]. 
The second term is 

(v tr (i))ik = + Kr(i))i* + ( 27 ) 

1 N 

1 i „ s 



( V tr(l))ik — Z • • • n • • • 5 iN k N 

x(i s |(D s (0)p s + p s D s (0))|A; s ), (28) 



v S(i))ik = ~TjY. 5 nki ■■■n--- 5 iN k N I du 5 + {uj) 

1V1 8=1 J -°° 

x (i s \([p s , v s }G s (uj + e is )~D s (uj) 

-T> s (lo)G s (lo + e is )[p s , v s ])\k s ) , (29) 

(int)\ 1 r s s ' r 

V tr(l))ik = 1^ d hki • • • l~l • • • l~l • • • d iN k N 

s<s' 

x ^(-l) P (^ s Pv|(D s (e Pis -e fc jp«' 
p 

+p s D s ,(ep i3l - e K ,)) \k s k s ,) , (30) 



1 N s f°° 

rr- Z • • • n • • • 5 iN k N J ^ du <5 4 



where 



D m (u) = -AixaZaiDi m (uj) , (31) 



ai (I = 1,2,3) are the Dirac matrices, Di m {uj) is the transverse part of the 
photon propagator in the Coulomb gauge. In the coordinate representation 
it is 

n /,x 1 ( exp(i\u\r) (exp(i\u\r) - 1) 1 

Afc(w,r) = - — ^ Oifc + ViVjfc = \. (32) 
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The contribution from the diagrams with two transverse and arbitrary num- 
ber of Coulomb photons is 

(U*r(2))*fc = Kr(2))ifc + ( u £(2))*k > ( 33 ) 

(i) i ^ 



2ttM s=1 

/oo 
gL>(? s |D s (cj)G s (cj + £i.)D s (a;)|fc s ) , (34) 
-oo 

Kr(2))jfc = 7^7- X! • • • l~l • • • l~l • • • 5 ijvfcjv 

s<s' 

x ]T(-i) p (P^Pv|D s (£p ls - £fc.) 

xD s ,(e Pis/ - e K/ )\k s k s ,) . (35) 

The formulas (23)-(35) were derived in [16]. The corresponding formulas 
for the case of a one-electron atom were first obtained in [14] (the overall sign 
of the contribution AE trty2 ) was corrected in [15,16]) and recently reproduced 
in [24,10]. 

The contributions v^ nt \ v^j^, and fj*™^ are leading for low aZ 
and completely define the nuclear recoil corrections within ^-(aZ) 4 approx- 
imation. It follows that within ^(aZ) 4 approximation the nuclear recoil 
corrections can be obtained by evaluating the expectation values with the 
Dirac wavefunctions of the operator 

R m = ^£(P S P*' -—{<*,+ r2 )Ps>) (36) 

s,s' s s 

In [25] the relativistic nuclear recoil corrections of order ^-(aZ) 4 to the 
energy levels of two- and three-electron multicharged ions were calculated 
using this operator. The expression (36) can be found by reformulating the 
Stone's theory as well [26]. 



3 Hydrogen-like atoms 

For hydrogen-like atoms the nuclear recoil corrections to the energy of a 
state a are defined by the diagonal matrix elements (AE = (v) aa ) of the 
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one-electrons contributions (24), (25), (28) ,(29), and (34). The terms AE*p 
and AE^y are leading at low Z. These terms can easily be calculated by 
using the virial relations for the Dirac equation [27-29]. Such a calculation 
gives [14] 



AE^ = 



to 



{» 



+ 



2(aZ) 4 



AE_ 



(i) 



tr(l) 

A£« 



2M l N 2 

+n r (4 7 2 -l)]}, 

to 2 {aZf 
~ If iV 4 7(4 7 2 - 1) 



AT4 7 (4 7 2 _ !) 



[k(2k(7 + n r ) - N) 



(37) 



[k(2k(7 + nr) - N) + n r (4 7 2 - 1)] , (38) 



TO 



(aZ) s 



2M 



M 2iV 2 



(39) 



where 



« = + 5) , 7 = V« 2 - (« Z ) 2 > 

iV = \Jn 2 — 2n r (|/t| — 7) , n = n r + \k\ , 

j is the total electron moment, I is the orbital moment, n is the principal 
quantum number, n r is the radial quantum number. Only these terms con- 

(«Z) 4 approximation. Expanding (39) in power series 



tribute within the ^ 
in aZ we find 



A£« 



to 2 c{aZf_ (aZy 
M 2n 2 + 2n 3 \? + 2 
(«Z) 6 a 



2n4( J + I)2V 4 ( J + i) 



(40) 



The terms AEf\ AE^ (1) , 



and 



the equations (25), (29), and 
(34) ) are given in the form that allows one to use the relativistic Coulomb 
Green function for their calculations. In addition, this form is convenient for 
aZ-expansion calculations [10]. However, in the present paper we transform 
these equations to ones that are most convenient for calculations using the 
finite basis set methods [30-32]. 

Integrating over uj in (25) we find 



A£f = -^£(a|p|n>Hp|a>. 

1VI e n <0 



(41) 
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(It should be noted here that the formula (41) was first found in [13]. Its 
derivation was refined in [14]. A similar formula but with the projector on 
the negative energy states of a free electron was obtained in the lowest order 
in aZ in [12].) The matrix elements of the momentum operator are easily 
calculated using the identity [25] 

p = - {aH + Ha) - aV c . (42) 
2 

Rotating in (29) the integration contour in the complex ui plane we find 
= AE%$ + AE™ + AES$ , (43) 
A4 2 (iJ = 2^ E {<a|p|n>(n|D(0)|a> 

+ (a|D(0)|n)(n|p|a)}, (44) 

AB« ( r p) = ~Ji Y. {{a\p\»)(»\U(e.-e n )\a) 

\Sn\<Sa 

+(a|D(e a -e n )|n)(ra|p|a)}, (46) 

S(y) = S 1 (y) + S 2 (y), 

G / x 7 exp (-yr) 
Si(y) = aZa , 

S 2 (y)=iaZ[Hj(y,r)n], 

exp (— yr)(l + yr) - 1 



where 



D(w) = Di(w)+D 2 (w), 

exp (i|tj|r) 
Di w = aZa , 

r 

D 2 (w) = iaZ[H, f(u, r)n] , 

1 — exp (i|a;|r)(l — i|a>|r) 



oo 2 r 2 
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_ , , ot iaZ rTT , 

D(0) = c*2- -— [if, n], 

n = ^. The term AE t ^j has real and imaginary parts. The imaginary part 
gives a small correction to the width of the level. Integrating over y in (45) 
and uniting the contributions AE t ^"j, AE t ^, and the real part of AE^j 
we find 

AE S(i) = ^f Re E (en-£a)(a|ia^(r)|n) 

x(n|[ia$i(r) + n$ 2 (r)]|a) , (47) 

where 

< = ^ + f , (48) 
*i(r) = — — [ci(A„r)sin(A n r) - si(A.„r) cos(A„r) 

'An'' 

/ x 71- ! , N 7r exp (iA n r) 

+sign(e a - e n - - 0(e a - e n Pl ; , 49 

2 A n r 

1 7T 

$ 2 (r) = -sign(e a - £„)— — — {-si(A n r) cos(A n r) - - + A n r 

(A n r)^ 2 

+ci(A„r)[sin(A n r) - (A„r) cos(A„r)] 

-(A n r) si(A„r) sin(A n r)} - ^ - #(e a - |e„|)7r/(A n , r) , (50) 

A n = |e a - £ n |, 9{x) = (x + |x|)/2x . 
The contribution AE tr ( 2 ) is equal 

AC = AiiS+^a + AEt 1 '' (51) 



A <>! - -^f*^ y » + X- e .)' <fl|s(,)|B)<B|s(,)|fl> - (52) 

AE^J = JL 2 <a|D(0)|n>(n|D(0)|a>, (53) 
A <?) = ^ E (a|D(e a -e B )|n>(n|D(e -e B )|a>- (54) 

|£n|<£a 

The term AE^j, like AE^j, has an imaginary part which gives a small 
contribution to the width of the level. 
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After integration over angles that is easily carried out using formulas 
presented in Appendix, the calculation of the expressions (41), (47), and 
(52)-(54) was done using the B-spline method for the Dirac equation, devel- 
oped in [31]. The zero boundary conditions and the grid selection algorithm 

4 

proposed in [33] were used. However, we used the grid r\ = where 

7o = \Jl — (aZ) 2 , instead of the grid r { = [33]. The radial integration 
caused no problems and was carried out with high accuracy using the Gauss- 
Legendre quadratures. The integration over y in (52) was also done by the 
Gauss- Legendre quadratures with a suitable transformation to map the in- 
finite integration range to a finite one. The uncertainty of the integration 
was estimated from the stability of the result with respect to change of the 
number of integration points and the grid parameters and was found to be 
much smaller than the uncertainty due to the finiteness of the basis set. The 
size of the box was chosen to be sufficiently large so as not to affect the 
results. The uncertainty, due to the finiteness of the basis set, was estimated 
by changing the number of splines from 40 to 90. In addition, to make an 
independent estimate of the uncertainty of the numerical results we calcu- 
lated the corrections AE^ and AE^L, using two different representations 

for them. So, the correction AE^ was calculated by the formula (41) as well 
as by 

AE& = -^{(a|p 2 |a) - £ (o|p|n>(n|p|a>}. (55) 

We found that the results of both calculations coincided with each other with 
good precision, and this coincidence improved when the number of splines 
increased. The correction AE^Ls was calculated by the equation (47) as well 
as by (43)-(45). The results of both calculations coincided with each other 
with high accuracy. 

Table 1, 2 and 3 show the results of the numerical calculation for the Is, 
2s, and 2pi states, respectively, expressed in terms of the function P(aZ) 
defined by 

AB» = A£f + AJ$„ + A4\> 2 , = ^P(»V (56) 

The functions P c , Pt r (i), and Pt r <2) correspond to the contributions AE^\ 
Ai^jL, and Ai^jL, respectively. For comparison, in the last columns of 
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the tables Salpeter's contributions [3-6] 



P$ s \aZ) 
P f s \aZ) 



-- In (aZ) - - 2.984129 + — In 2 + — 
3 V ; 3 3 9 



(57) 



(2pi) 




(58) 



(59) 



are given. The uncertainties given in the tables correspond only to errors 
of the numerical calculation. In addition, there is an uncertainty due to 
deviation from the point single particle model of the nucleus, used here. 

To make a more detailed comparison with the aZ-expansion calculations 
we represent the functions P c , Ptr(i), and Pt r (2) for the s states in the form 

P c = ai + a 2 aZ + a 3 (aZ) 2 ln(aZ) + a A (aZ) 2 , 
-Ptr(i) = h In (aZ) + b 2 + b 3 aZ In (aZ) 

+b A aZ + b 5 (aZ) 2 In (aZ) + b 6 {aZ) 2 + b 7 (aZ) 3 , 
Ptr{2) = ci In (aZ) + c 2 + c 3 aZ In (aZ) 

+c 4 aZ + c 5 {aZ) 2 In (aZ) + c 6 (aZ) 2 + c 7 («Z) 3 . (60) 

The coefficients a^, 6j, and q can be calculated from our numerical results 
for the P(aZ)-functions. Such a calculation for the 2s state using the values 
of the P ( a Z) -functions for Z =1,2,3,5,8,15,30 gives 

a\ = —1.3333, a2 = 3.156, 

6i = -2.6662, 6a = -0.091, 6 3 = -6.02, 6 4 = -9.98, 

ci = 2.0031, c 2 = 4.338, c 3 = 6.46 , c 4 = 5.92. (61) 

The coefficients ai, 6i ;2 , and Ci )2 are in good agreement with Salpeter's results 



Within errors of the numerical procedure our values 63, C3, are in good agree- 
ment with the analytical result of [8,9] 



ai = -1.3333 61 = -2.6666, b 2 = -0.094, 
ci = 2.0000, c 2 = 4.318. 



(62) 



c 3 = -2tt = -6.2832 , b 3 + c 3 = 



(63) 
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(the coefficient 63 was first found in [7]). The coefficient a 2 coincides, within 
the numerical errors, with the corresponding coefficient (02 = 7r = 3.1459) 
obtained in [7]. The coefficients 64 and c 4 are in satisfactory agreement with 
the results of [10] 

b A = -10.996, c 4 = 5.569. (64) 

For the Is state we have found a similar agreement. 

To make a similar comparison for the 2p± state we represent the functions 
P tr (i) and P tr (2) for this state in the form 

Ptr(i) = h + b 2 aZ + b 3 (aZ) 2 \n(aZ) 

+h(aZ) 2 + b 5 (aZ) 3 In (aZ) + b 6 (aZ) 3 + b 7 (aZ) A , 
Ptr(2) = ci + c 2 aZ + c 3 (aZ) 2 In (aZ) 

+c 4 (aZ) 2 + c 5 (aZf In (aZ) + c Q (aZf + c 7 (aZ) 4 . (65) 

Using our values of P(aZ) for Z — 1, 2, 3, 5, 8, 15, 30 we have found 

h = -0.142178, 6 2 = -0.26166, 

Cl = -0.166666, c 2 = 1.30881. (66) 

The coefficients bi and c x are in excellent agreement with the Salpeter's 
results: 61 = -0.142178 and c x = -0.166667. Adding to the sum b 2 + c 2 
the corresponding coefficient from the equation (40), we find that the total 
coefficient of the ^ ^"f^ contribution for the 2pi state is 1.43985. The 

related analytical result obtained in [11] is |j7r = 1.43990. 

The term AE^ does not contribute to the Lamb shift of hydrogen-like 
atoms. The contribution of the difference between AE^ and the Salpeter's 
correction to the Lamb shift (n = 2) of hydrogen is —1.32(6) kHz. The 
corresponding result for the ground state is —7.1(9) kHz. These results are 
in good agreement with analytical calculations of the ^-(«Z) 6 contributions 
[10,11]. So, according to [10] the total ^(aZ) 6 correction, including the 
related term from the equation (40), is —7. 4kHz and —0.77kHz for the Is 
and 2s states, respectively. The ^-(aZ) 6 correction for the 2j9i state, found 

in [11], is 0.58 kHz. (We note that in [11] the correction of order ^(a) 2 («Z) 4 
for p states is also calculated.) 
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Let us consider the nuclear recoil corrections for hydrogen-like uranium. 
According to the formula (39) the first correction is 

AE$ = 0.26 eV , AE$ = AE^ = 0.08 eV . (67) 

The second correction defined by (56) is 

AE { $ = 0.24 eV , AE^f = 0.05 eV , AE^ = 0.01 eV . (68) 

2 

In the next section we use these results to find the total nuclear recoil con- 
tribution to the energy of the 2pi — 2s transition in lithium-like uranium. 



4 High Z lithium-like atoms 

The wavefunction of a high Z lithium-like atom with one electron over the 
closed (Is) 2 shell in the zeroth approximation is 

u = ^=f E(- 1 ) P ^t(^ P1 )V'i s |( J P2)^(P3) . (69) 

The nuclear recoil correction for the lithium-like atom is the sum of the 
one- and two-electron corrections. The one-electron correction is obtained 
by summing all the one-electron contributions considered in the preceeding 
section over all the one-electron states that are occupied. According to (26), 
(30), and (35) the two-electron corrections for the state considered here are 



AE™ = ~ E (a\p\n)(n\p\a) , (70) 
j-f E {(a|p|n)(ra|D(e -e n )|a) (71) 

£n— £ls 

+ (a\H(e a - e n )\n)(n\p\a)} , 



A rr( int ) . 



^$5) = ~ E (a|D( eo - en )|n>(n|D( ea - en )|a>. (72) 



1 



£n=£l 



The terms AE^f} and AE^l have real and imagine parts and are cancelled 
by a part of the one-electron terms AE[^ and AE^'^), which corresponds 
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to the Is states. So, for the (ls) 2 2s and (ls) 2 2pi states the imagine parts of 
the one- and two-electron contributions are completely cancelled. 

We note here that the nuclear recoil corrections for a high Z lithium-like 
atom with one electron over the closed (Is) 2 shell can be obtained from the 
nuclear recoil corrections for the hydrogen-like atom by changing the sign 
of iO in the denominators of the electron propagator in the Coulomb field 
of the nucleus, corresponding to the states of the closed shell. It follows, in 
particular, the sum of the one- and two-electron Coulomb contributions can 
be represented in a simple form 

AE < = v^{ E IHpNI 2 - E IHpHI 2 }- (73) 

e n >sis e n <s la 

The table 4 shows the results of the calculation of the corrections (70), (71), 
and (72) for the (ls) 2 2pi state (for the (ls) 2 2s states these corrections are 
equal to zero), expressed in terms of the function Q(aZ) defined by 

AE int = AE (int) + + A £M) = _^ {aZ) 2 Q{aZ) (?4) 

Here we have taken into account the known non-relativistic limit of this 
correction [34]. Within the ^-(aZ) 4 approximation the function Q[aZ) that 
we denote by Ql{&Z) is [25] 

Q L {aZ) = l + (aZ) 2 (-- + ln-). (75) 

For comparison, this function is given in the table as well. The functions 
Q c (aZ), Qtr^i^aZ), and Q tr (2)(o:Z) correspond to the corrections AE^ nt \ 

AE^f}, and AE^%}, respectively. In leading orders in aZ they are 



Q c {aZ) = l + («Z) 2 (| + ln|), (76) 



Q tr{1) (aZ) = - 7 -(aZ) 2 , (77) 
49 

Q tr[2) (aZ) = ^{aZf. (78) 

For low Z, in addition to the corrections considered here, the Coulomb 
electron-electron interaction corrections to the non-relativistic nuclear recoil 
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contribution must be calculated separately. The main contribution from 
these corrections is of order ^(aZ) 2 ^. 

Sometimes, to estimate the nuclear recoil corrections for high Z the non- 
relativistic nuclear recoil operator is averaged with the Dirac wavefunctions. 
But, as one can see from the formulas (75)-(77) and the table 4, like the one 
electron case (see the formulas (37)-(40)), this contribution is considerably 
cancelled by the one-transverse-photon contribution. 

According to [35] the experimental value of the energy of the (ls) 2 2pi — 
(ls) 2 2s transition in lithium-like uranium is 280.59(10) eV^. Let us find the 
total nuclear recoil contribution to the energy of this transition. According 
to our calculation the term AE mt is — 0.03 eV. Adding to this value the 
one-electron contribution defined by (68) we find 

A£ ( i s) 2 2pi - AE {ls) 2 2s = -0.07 eV. 

This correction, largely made up of the QED contributions, is comparable 
with the uncertainty of the experimental value and, hence, will be important 
for comparison of theory with experiment, when calculations of all diagrams 
in the second order in a are completed. 
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Appendix 

The integration over angles in the expressions considered here is carried out 
using the formula 

^(n 1 j 1 l 1 m 1 \A\n 2 j2l2rn2)(n2j2l2rn 2 \B\n 1 j 1 l 1 m 1 ) 

1712 

= { _ l)h+j2 -2m, ^^(n^Z^l^llnajaZaJCnajaZall^llnjjiZi) , (79) 

where (riiji^i| \A X \ \ri2j2h) , (^2.72^2! I-B 1 ! I^iji^i) are the reduced matrix ele- 
ments [36]. For A = a:0(r) , n0(r) one can find 



-(-1)% 2 { j } jf /n 1Jlil (r)^ t2 , 2i2 (r)0(r)r 2 rfr] , (80) 

1 r ■ • f°° 

(n 1 j 1 Z 1 ||n0(r)||n 2 j 2 / 2 ) = (-1) J2 ~2 [Z/^ 2 jf c/„ ljl i 1 (r)^ n2j2Z2 (r)0(r)r 2 dr 

+Zgg fn inh (r)fn 2j2 i 2 (r)<P(r)r 2 dr], (81) 



where 



^V 2 2 = V( 2/ i + !)(2'2 + l)(2ji + l)(2j 2 + 1) qJ I^V^] , (82) 

/' = 2j — /; g n jiij-) and f n ji( r ) are the upper and lower radial components of 
the Dirac wavef unction [37]: 

( r \ _ ( 9nji(r)£l jlm (n) \ 
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Table 1: The results of the numerical calculation of the one-electron nuclear 
recoil corrections to the Is state energy expressed in terms of the function 
P(aZ) defined by equation (56). Ps{olZ) is the Salpeter's contribution de- 
fined by equation (57). 



z 


P c (aZ) 


P tr{ i)(aZ) 


P tr(2 )(aZ) 


P(aZ) 


P s (aZ) 


1 


-1.3111(2) 


12.568(2) 


-5.8267(3) 


5.430(2) 


5.4461 


5 


-1.2345(1) 


8.5854(3) 


-3.0476(2) 


4.3033(4) 


4.3731 


10 


-1.1586 


6.9974(1) 


-2.0438 


3.7950(1) 


3.9110 


15 


-1.0994 


6.1340(1) 


-1.5373 


3.4973(1) 


3.6407 


20 


-1.0537 


5.5678(1) 


-1.2201 


3.2940(1) 


3.4489 


25 


-1.0192 


5.1671(1) 


-0.9996 


3.1483(1) 


3.3001 


30 


-0.9946 


4.8744(1) 


-0.8362 


3.0437(1) 


3.1786 


35 


-0.9790 


4.6598(1) 


-0.7094 


2.9714(1) 


3.0758 


40 


-0.9721 


4.5065(1) 


-0.6076 


2.9268(1) 


2.9868 


45 


-0.9740 


4.4048(1) 


-0.5231 


2.9077(1) 


2.9083 


50 


-0.9849 


4.3496(1) 


-0.4510 


2.9137(1) 


2.8380 


55 


-1.0059 


4.3389(1) 


-0.3874 


2.9456(1) 


2.7745 


60 


-1.0383 


4.3739(2) 


-0.3295 


3.0061(2) 


2.7165 


65 


-1.0845(1) 


4.4588(2) 


-0.2746 


3.0997(2) 


2.6631 


70 


-1.1479(2) 


4.6014(3) 


-0.2201 


3.2334(4) 


2.6137 


75 


-1.2339(3) 


4.8153(7) 


-0.1631 


3.4183(8) 


2.5677 


80 


-1.3506(5) 


5.122(1) 


-0.0996(1) 


3.672(1) 


2.5247 


85 


-1.512(1) 


5.558(4) 


-0.0237(2) 


4.022(4) 


2.4843 


90 


-1.741(3) 


6.186(7) 


0.0743(9) 


4.519(8) 


2.4462 


92 


-1.861(5) 


6.51(1) 


0.123(1) 


4.77(1) 


2.4315 


95 


-2.084(9) 


7.12(3) 


0.212(1) 


5.25(3) 


2.4101 


100 


-2.64(3) 


8.6(1) 


0.428(6) 


6.4(1) 


2.3759 
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Table 2: The results of the numerical calculation of the one-electron nuclear 
recoil corrections to the 2s state energy expressed in terms of the function 
P(aZ) defined by equation (56). Ps{olZ) is the Salpeter's contribution de- 
fined by equation (58). 



z 


P c (aZ) 


Ptr(i)(aZ) 


P tr{ 2)(aZ) 


P(aZ) 


P s (aZ) 


1 


-1.3112(2) 


13.177(1) 


-5.7103(3) 


6.155(1) 


6.1710 


5 


-1.2351(1) 


9.1911(2) 


-2.9225(1) 


5.0335(2) 


5.0980 


10 


-1.1612 


7.6075(1) 


-1.9080 


4.5383(1) 


4.6359 


15 


-1.1055 


6.7562 


-1.3908 


4.2599 


4.3656 


20 


-1.0647 


6.2093 


-1.0621 


4.0825 


4.1738 


25 


-1.0367 


5.8352 


-0.8294 


3.9691 


4.0251 


30 


-1.0202 


5.5767 


-0.6528 


3.9037 


3.9035 


35 


-1.0147 


5.4047 


-0.5115 


3.8785 


3.8008 


40 


-1.0202 


5.3037 


-0.3935 


3.8900 


3.7117 


45 


-1.0372 


5.2656 


-0.2908 


3.9376 


3.6332 


50 


-1.0668 


5.2876(1) 


-0.1980 


4.0228(1) 


3.5630 


55 


-1.1108 


5.3711(1) 


-0.1105 


4.1498(1) 


3.4994 


60 


-1.1723(1) 


5.5218(1) 


-0.0247 


4.3248(2) 


3.4414 


65 


-1.2554(1) 


5.7504(2) 


0.0634 


4.5584(2) 


3.3881 


70 


-1.3668(2) 


6.0743(4) 


0.1581(1) 


4.8656(5) 


3.3387 


75 


-1.5164(4) 


6.5211(7) 


0.2651(1) 


5.2698(8) 


3.2927 


80 


-1.7199(7) 


7.135(2) 


0.3921(2) 


5.807(2) 


3.2496 


85 


-2.003(1) 


7.988(4) 


0.5516(4) 


6.537(4) 


3.2092 


90 


-2.413(4) 


9.205(8) 


0.7645(6) 


7.557(9) 


3.1711 


92 


-2.630(7) 


9.84(1) 


0.872(1) 


8.08(2) 


3.1565 


95 


-3.04(2) 


11.02(2) 


1.070(2) 


9.05(3) 


3.1351 


100 


-4.07(5) 


13.9(1) 


1.55(1) 


11.4(2) 


3.1009 
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Table 3: The results of the numerical calculation of the one-electron nuclear 
recoil corrections to the 2pi state energy expressed in terms of the func- 
tion P(aZ) defined by equation (56). Ps(aZ) is the Salpeter's contribution 
defined by equation (59). 



z 


Pc(aZ) 


Ptr(i)(aZ) 


P tr(2) (aZ) 


P(aZ) 


P s {ctZ) 


1 


-0.0000 


-0.1440 


-0.1571 


-0.3011 


-0.3088 


5 


-0.0007 


-0.1492 


-0.1194 


-0.2692 


-0.3088 


10 


-0.0024 


-0.1526 


-0.0727 


-0.2277 


-0.3088 


15 


-0.0051 


-0.1535 


-0.0258 


-0.1845 


-0.3088 


20 


-0.0088 


-0.1524 


0.0218 


-0.1393 


-0.3088 


25 


-0.0133 


-0.1493 


0.0706 


-0.0920 


-0.3088 


30 


-0.0189 


-0.1444 


0.1212 


-0.0421 


-0.3088 


35 


-0.0255 


-0.1375 


0.1742 


0.0112 


-0.3088 


40 


-0.0335 


-0.1284 


0.2304 


0.0685 


-0.3088 


45 


-0.0432 


-0.1165 


0.2906 


0.1310 


-0.3088 


50 


-0.0548 


-0.1012 


0.3560 


0.2000 


-0.3088 


55 


-0.0691 


-0.0814 


0.4278 


0.2774 


-0.3088 


60 


-0.0868 


-0.0555 


0.5078 


0.3655 


-0.3088 


65 


-0.1091 


-0.0211 


0.5982 


0.4680 


-0.3088 


70 


-0.1376 


0.0252 


0.7018 


0.5894 


-0.3088 


75 


-0.1750 


0.0891(1) 


0.8229 


0.7370(1) 


-0.3088 


80 


-0.2253(1) 


0.1796(1) 


0.9671 


0.9214(2) 


-0.3088 


85 


-0.2954(2) 


0.3123(3) 


1.1429(1) 


1.1598(4) 


-0.3088 


90 


-0.3972(6) 


0.515(1) 


1.3632(1) 


1.481(1) 


-0.3088 


92 


-0.451(1) 


0.626(1) 


1.468(2) 


1.643(3) 


-0.3088 


95 


-0.554(2) 


0.842(3) 


1.649(3) 


1.937(5) 


-0.3088 


100 


-0.816(9) 


1.41(1) 


2.040(3) 


2.63(2) 


-0.3088 
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Table 4: The results of the numerical calculation of the two-electron nuclear 
recoil corrections AE^ int ^ for the (ls) 2 2pi state of lithium-like ions expressed 
in terms of the function Q(aZ) defined by equation (74). Q L (aZ) is the 
leading contribution defined by equation (75). 



z 


QJaZ) 


Qtrl^id-Z) 
LI \i) \ / 


QtrC2)(o:Z) 

LI \6) \ J 


Q(aZ) 


QdaZ) 


5 


1.00168 


-0.00233 


0.00000 


0.99935 


0.99935 


10 


1.00677 


-0.00938 


0.00002 


0.99741 


0.99741 


15 


1.01533 


-0.02129 


0.00011 


0.99416 


0.99417 


20 


1.02753 


-0.03830 


0.00036 


0.98959 


0.98964 


25 


1.04359 


-0.06077 


0.00088 


0.98370 


0.98381 


30 


1.06378 


-0.08920 


0.00186 


0.97645 


0.97669 


35 


1.08851 


-0.12422 


0.00353 


0.96782 


0.96827 


40 


1.11827 


-0.16669 


0.00617 


0.95776 


0.95856 


45 


1.15370 


-0.21767 


0.01019 


0.94622 


0.94755 


50 


1.19560 


-0.27853 


0.01607 


0.93313 


0.93525 


55 


1.24500 


-0.35105 


0.02447 


0.91841 


0.92165 


60 


1.30322 


-0.43751 


0.03625 


0.90195 


0.90676 


65 


1.37198 


-0.54091 


0.05254 


0.88361 


0.89057 


70 


1.45352 


-0.66521 


0.07488 


0.86320 


0.87309 


75 


1.55087 


-0.81573 


0.10538 


0.84052 


0.85431 


80 


1.66810 


-0.99980 


0.14699 


0.81529 


0.83424 


85 


1.81092 


-1.22771 


0.20395 


0.78716 


0.81287 


90 


1.98751 


-1.51431 


0.28250 


0.75570 


0.79021 


92 


2.07014 


-1.65003 


0.32196 


0.74206 


0.78078 


95 


2.21001 


-1.88186 


0.39221 


0.72035 


0.76625 


100 


2.49719 


-2.36503 


0.54826 


0.68041 


0.74099 
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